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ABSTRACT 

The  research  performed  under  AFOSR  Grant  87-007^  during  the  first  year 
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flow  control  via  an  unsteady  free  stream. 
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1 .  RESEARCH  OBJECTIVES 

The  objectives  that  were  set  for  this  grant  consisted  of  the  following: 
(i)  To  examine  and  achieve  active  control  of  two-dimensional  (2-D) 

unsteady  separated  flows  through  direct  numerical  simulation  using 
the  unsteady  Navi er-Stokes  (MS)  equations, 


(ii)  To  study  three-dimensional  (3-D)  unsteady,  separated  vortex- 

dominated  flows  by  developing  a  direct  implicit  technique;  this 
should  also  permit  analysis  of  flow  instabilities  using  model 


problems . 

The  research  pursued  in  both  of  these  areas  since  October  1986  and  the 
progress  achieved  is  described  briefly  in  the  next  section. 


2.  RESEARCH  PROGRESS 


2 . 1  Control  of  Separated  Flows 

Active  control/management  of  separated  flows  is  studied  by  considering 
the  flow  in  a  backstep  channel  with  an  oscillating  flap  located  near  the 
downstream  end  of  the  primary  separation  bubble  on  the  lower  wall  as  shown 
in  Fig.  1.  Time-dependent  boundary-or iented  clustered  conformal  coordinates 
were  generated,  during  the  first  half  of  this  reporting  period,  using  the 
Schwarz-Christof fel  conformal  mapping,  followed  by  suitable  contraction 
mappings  to  resolve  the  critical  regions  in  this  flow.  A  brief  description 
of  the  method  for  developing  the  transformation  and  determining  the 
coordinates  for  arbitrary  orientations  of  the  oscillating  flap  is  given  in 
Appendix.  During  the  second  half  of  this  reporting  period,  effort  has  been 
focussed  on  determination  of  the  flow  in  the  temporally  deforming  domain  of 


this  flow. 
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2.1.1  Governing  Equations 


f  ur.c  t  i  o  n  formulation 


Incompressible  flow 


‘‘rnros^iitsc  inst  hsnat  i  os  1 1  v  by  bh*?  vorti  city-  stream 
of  the  time-dependent  Gavier-GtCr.es  equations  for  1-1 


equations  for 


unsteady  flow,  the  stream  function  v  is  obtained  oy  a  direct-solution  method 
based  on  bloek-Gaussian  elimination  'SGE> ,  for  every  time-step  advance  of 
the  vorti city  field  u.  With  temporally  deforming  grids,  even  the  linear 
differential  equation  for  y  takes  on  time-varying  coeff  icie.nts .  Therefore, 
additional  computational  time  is  required  not  only  for  the  evaluation  of 
these  coefficients  but  also  for  the  solution  of  the  corresponding  partial 
differential  equation  by  direct  inversion.  3y  contrast,  for  efficiency  of 
the  3GE  procedure,  the  differential  operator  for  tj;  must  remain  time- 
invariant.  Conformal  mapping  allows  this  to  be  achieved  even  in  the 
presence  of  temporally  deforming  grids,  as  briefly  described  next. 

In  Cartesian  coordinates,  the  stream-function  equation  has  the  form 

(1) 

xx  yy 

where  the  source  function  is  the  vorticity  <u  which  is  a  function  of  x,  y  and 
t ,  i .e  .  , 

oj  =  u>(x  ,y  ,t )  .  (2) 

In  conformal  coordinates  (£, n) »  this  equation  can  be  represented  as 


where 


*5£  *  *r,n  ■  •  h  “ 


h  =  h  =  h. 


is  the  scale  factor  of  the  conformal  mapping 
£  +  in  =  f ( x- iy ) 


In  terms  of  complex  variables,  the  conformal  mapping  and  its  scale  factor 


become 

z  =  fU) 

where 


h 


Jdf 

dz 

|d? 

z,  =  5  +  in  and  z  =  x  +  iy 


(6a , b ) 


(7a, b) 


For  a  temporally  deforming  boundary-al i gned  conformal  coordinate 
system,  the  scale  factor  h  varies  with  time,  in  addition  to  being  spatially 
varying.  However,  by  virtue  of  the  special  form  of  the  stream  function 
equation  (3)  that  results  with  a  conformal  mapping,  this  time-varying  h  can 
be  absorbed  in  the  source  function  of  this  equation  leading  simply  to  the 
equation 


where 


^  +  %n  =  GU’n’t: 


GU,n,t)  =  -  h2cd  . 


(3a) 


(8b) 


This  leaves  the  stream-f unction  equation  with  time-invariant  coefficients. 

It  is  important  to  note  that  any  clustering  transformations,  introduced  in 
Eq .  (8a)  for  the  purpose  of  resolving  the  length  scales  of  the  problem,  must 
be  time- independent  in  order  for  the  final  form  of  the  equation  to  possess 
time- independent  coefficients.  The  clustering  transformations  used  in  the 
present  work  are  of  the  form  £  =  £(S)  ,  n  =  n(N)  where  (S,N)  are  the  final 
computational  coordinates. 

The  transport  equation  for  the  vorticity  u  far  2-D  incompressible  flow 


is  written  in  vector  form  as 


a).  +  V*Vto  =  —  V  a) 
t  Re 

Here,  V  is  the  velocity  vector  with  components  u  and  v  along 

coordinate  directions  E,  and  n,  respectively,  related  to  the  s 

<Jj  as  follows: 

1_  _  1  3  vp 

U  ~  h  3n  :  V'~h3l  * 

The  time-dependent  conformal  mapping  can  be  expressed  as 
E  =  E(x,y,t) 
n  =  n(x,y  ,t ) 

T  =  t 


'  9 


the  cor.formal- 
tream  functioi 


(10) 


where  t  is  the  tranformed  time  coordinate.  Hence,  the  time  derivative  tern 
uit  in  Eq.  (9)  transforms  as 


'JJ.  =  0)  +  ;,D.  +  h  'JO 

t  t  t  (  t  n 


1  ) 


or,  finally,  as 


1 


=  a)  +  — 5  )(x -  y_x  )ojc  +  (xry  -  y,x  )w„]  . 
h‘ 


t  ,  2  — nJ  t  7  J  /wrr 


(12) 


The  spatial  derivatives  of  x,  y  are  easily  evaluated  using  the  complex 


quantity  dz/d £.  The  grid  speed  terms  x  ,  y ^  are  determined  using  the 


mapping  at  two  successive  levels  of  time. 

2.1.2  Boundary  Conditions 

The  wall  boundary  conditions  are  obtained  from  the  condition  of  zero 
slip  along  the  non-porous  walls  of  the  backstep  channel.  These  are  then 
expressed  in  terms  of  the  stream  function  ■}>,  as  follows: 

Along  n  =  n  .  and 

mi  n 


1 

i 

1 

i 


n 


For  points  on  a  stationary  non-porous  wall,  these  lead  to  the  conditions 


that,  at  these  points, 

=0  and  |  =  0  ( 1 4a  ,b) 

k  nk 

where  ns  denotes  stationary-boundary  points.  Equation  (14a)  implies  that 

the  boundary  is  a  streamline.  For  points  on  a  moving  boundary,  the  right- 
hand  sides  of  either  Eq.  (13a)  or  Eq.  (13b)  or  both  may  be  non- zero,  even 
when  the  boundary  is  non-porous.  This  is  easily  confirmed  by  examining  the 
streamline  pattern  around  a  cylinder  moving  at  a  constant  velocity  in  an 
otherwise  stagnant  fluid  (Fig.  2).  If,  at  points  A  and  9  on  the  cylinder, 

at  A  and  ^  ^  at  B,  then  ^  is  not  constant  on  the  cylinder  surface. 

Hence,  for  the  present  problem,  >4  along  the  moving  flap  is  no  longer  zero 
even  though  it  is  zero  elsewhere  on  the  lower  boundary  of  the  channel.  The 
distribution  of  \\>  along  the  moving  flap  is  obtained  by  using  Eq .  (13) 
together  with  information  about  the  flap  motion,  i.e., 

Along  the  flap: 

u„  =  0  and  v  =  re  ( 1 5a  ,b ) 

r  r 

where  r  is  the  distance  along  the  flap  measured  from  its  base  (Fig.  3).  and 

9  is  the  instantaneous  rate  of  its  rotation.  The  specified  motion  of  the 

flap  will  provide  9,  whereas  r  is  obtained  from  £  and  h  f or  the  conformal 
mapping.  Use  of  Eqs.  (15a,b)  in  Eqs.  ( 1 3a , b )  leads  to 

^2 

4 u  -  =  -  S  h  v  d£  and  4j  I  =0  (I6a,b) 

2  1  _  F  n  n- 

*1  F 

As  regards  concerns  about  accuracy  of  the  integration  indicated  in  Eq. 

( 1 6a ) ,  it  is  important  to  note  that  41  =  0  at  points  adjoining  both  sides  of 
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the  base  of  the  flap.  This  provides  a  check  on  the  integration  accuracy  and 
is  used  to  formulate  the  discretized  integration. 

Another  concern  that  needs  to  be  dealt  with  is  the  following.  The  grid 
points  in  the  computational  plane  are  uniformly  distributed  and  fixed  for 
all  time.  With  time-dependent  clustering  transformations,  the  grid-point 
distribution  in  the  doubly-i nf i nite  straight-channel  plane  is  also  fixed  for 
all  time.  The  corresponding  distributions  along  the  flap  in  the  physical 
plane  at  two  successive  time  levels  may  not,  in  general,  be  related  by  the 
physical  motion  of  the  flap;  see  Fig.  4.  For  a  straight  flap,  this  can  be 
treated  without  difficulty,  using  simple  interpolation.  For  curved  flap 
resulting  from  possible  binding,  this  requires  careful  consideration. 

All  of  the  above  analysis  has  been  incorporated  into  the  flow  program 
and  is  being  made  operational.  Some  preliminary  results  are  anticipated 
within  the  next  few  weeks,  unless  any  additional  problems  are  encountered 
and  need  to  be  resolved.  Current  plans  are  for  a  paper  based  on  this  work 
to  be  presented  at  the  National  Fluid  Dynamics  Congress  scheduled  for  July, 
1983,  in  Cincinnati,  Ohio, 

2.1.3  Flow  Control  for  Axisymmetric  Configurations 

For  2-D  rectangular  configurations,  the  stream  function  equation 
possesses  constant  coefficients.  This  is  not  so  for  axisymmetric 
configurations  where  some  of  the  coefficients  in  the  equation  governing  the 
Stokes  stream  function  contain  the  radius  r  in  an  essential  manner.  For  a 
given  computational  point,  the  value  of  r  changes  with  time  as  the  boundary 
shape  alters.  This  implies  that  the  Schwar z-Chr i stof f el  mapping  obtained 
directly  as  dz/dc  needs  to  be  integrated  at  every  time  step.  Furthermore, 
this  implies  that  the  elimination  phase  of  the  8GE  method  needs  to  be 
repeated  at  every  time  step.  The  time-consuming  nature  of  this  phase  has  a 
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large  adverse  effect  on  the  efficiency  of  the  overall  procedure.  The 
possible  use  of  the  MG-SI  (Multi-Grid  Strongly  Implicit)  procedure  developed 
earlier  by  U.  Ghia  et  al .  [1983.  1987]  is  being  explored  for  the  stream 
function  equation  for  the  axisymmetric  configuration  of  a  pipe-orifice  with 
bending  of  the  orifice  plate.  The  preliminary  findings  from  this  effort  are 
summarized  in  Table  I,  and  indicate  that  the  MG-SI  scheme  is  a  viable 
procedure  for  solution  of  elliptic  equations  with  time-dependent 
coefficients . 

2.2  Unsteady  Free  Stream  -  A  Control  Mechanism 

Aerodynamic  forces  can  be  significantly  affected  by  the  unsteady  nature 
of  the  separation  and  reattachment  phenomena  on  the  suction  surface  of  an 
aerodynamic  body  at  high  angle  of  attack.  Present  understanding  of  the 
underlying  mechanisms  is  quite  limited.  A  study  of  the  flow  past  an 
elliptic  cylinder  at  angle  of  attack  was  initiated  around  July  1987,  with 
the  objective  of  studying  flow  control  by  subjecting  the  elliptic  cylinder 
to  a  time-dependent  free  stream. 

As  a  first  phase  of  this  work,  the  flow  past  an  elliptic  cylinder  is 
being  studied,  using  the  two-dimensional  unsteady  incompressible  Navier- 
Stokes  analysis  of  K.  Ghia,  Osswald  and  U.  Ghia  [1985].  Clustered  conformal 
transformations  are  developed  to  obtain  appropriate  surface-oriented 
coordinates  for  an  elliptic  cylinder.  The  grid-clustering  transformations 
implemented  attempt  to  resolve  the  multiple  disparate  length  scales  of  the 
unsteady  separated  flow;  a  typical  grid  consists  of  (229x45)  mesh  points. 

The  results  of  the  present  unsteady  analysis  are  compared, 
qualitatively,  with  the  experimental  data  of  Izumi  and  Kuwahara  [1983]  in 
Figs.  5  and  6.  For  the  case  of  an  elliptic  cylinder,  with  50%  thickness, 


angle  of  attack  a  =  20°  and  Reynolds  number  Re  =  369.  Before  discussing 
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these  figures,  it  is  appropriate  to  point  out  the  conditions  under  which  the 
two  sets  of  data  are  obtained,  so  the  qualitative  comparison  carried  out  is 
in  proper  perspective.  In  the  experiments,  an  electrolytic  precipation 
method  was  used  to  obtain  flow  visualization  results;  Taneda  [1977]  has 
shown  that  such  an  approach  leads  to  visualization  of  'integrated 
streaklines'  which,  in  turn,  show  the  'vorticity'  layers.  Further,  although 
the  time  scales  are  non-dimensional ized  in  the  same  way,  the  cylinder  was 
accelerated  over  a  finite  amount  of  time  in  the  experiment,  whereas  in  the 
present  analysis,  the  cylinder  was  impulsively  started  from  rest.  Finally, 
the  orientation  of  the  two  sets  of  results  shown  is  also  different;  the 
experimental  photographs  show  a  horizontal  free  stream  and  the  cylinder  at 
an  angle  of  attack,  but  the  present  computational  results  are  plotted  with 
the  axis  of  the  cylinder  being  horizontal  and  the  oncoming  flow  is  at  the 
appropriate  angle  of  attack. 


In  Fig.  5,  the  instantaneous  flow  field  results  for  "vorticity"  as 
stated  above,  are  shown  for  characteristic  time  T  =  0.86.  In  the  photograph 
as  well  as  the  plot,  a  counterclockwi se  rotating  vortex  is  beginning  to  shed 
off  of  the  lower  surface  near  the  trailing  edge,  while  a  clockwise  rotating 
vortex  sheet  is  reaching  across  the  upper  surface  and  curls  down  towards  the 
trailing  edge.  In  Fig.  6,  the  results  are  depicted  at  T  =  1.69.  Both  sets 
of  results  show  that  the  counterclockwise  rotating  vortex,  that  was 
initially  seen  being  shed  at  T  =  0.86,  has  moved  downstream  away  from  the 
cylinder.  Figure  6  also  shows  the  vortical  structure  positioned  just  over 
the  upper  portion  of  the  trailing  edge  and  the  vortex  sheet  from  the  leading 
edge  curling  up  and  over  this  vortical  structure  in  both  the  photograph  as 
well  as  the  plot  of  present  results. 
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Preliminary  results  are  obtained  for  flow  past  an  elliptic  cylinder 


with  25%  thickness  at  an  angle  of  attack  a  =  15°  and  Re  =  1000.  The  flow 
solutions  monitored  showed  that  the  cylinder,  which  is  impulsively  started 
from  rest,  goes  through  a  highly  unsteady  phase  until  T  =  17,  after  which 
the  asymptotic  state  in  the  form  of  a  limit  cycle  is  reached.  Further 
detailed  examination  of  the  limit  cycle  solution  showed  alternate  shedding 
of  vortices  from  the  suction  and  pressure  surfaces.  The  instantaneous 
stream  function  and  vorticity  are  plotted  in  Fig.  7  (a-d)  for  the  above 
configuration  at  T  =  28.47  corresponding  to  the  maximum-lift  instant  and 
T  =  29.106  for  the  minimum-lift  instant  in  the  same  cycle.  The 
corresponding  instantaneous  velocity  vectors  for  the  two  cases  are  depicted 
in  Fig.  7  (e-f ) . 

2. 3  Three-Dimensional  Incompressible  Separated  Flow 

This  study  has  continued  using  the  model  problem  of  the  shear-driven 
flow  within  a  3-D  rectangular  cavity.  Preliminary  results  using  the  direct 
inversion  methodology  for  the  unsteady  Naver-Stokes  equations  were  obtained 
during  the  first  half  of  this  report  period  culminating  in  the  paper 
presented  by  Osswald,  K.  Ghia  and  U.  Ghia  [1987]  at  the  8th  AIAA  CFD 
Conference,  Honolulu,  Hawaii,  appearing  in  the  conference  proceedings. 
Results  obtained  on  the  Department's  Perkin  Elmer  3250  using  the  serial 
research  version  of  the  code  were  presented  for  a  (17x17x17)  grid.  This 
proved  the  formulation  of  the  unsteady  Navier-Stokes  problem  as  well  as  the 
efficiency  of  the  direct,  inversion  numerical  procedures. 

Since  then  the  code  has  been  fully  vectorized  for  the  CYBER  205 
computer,  resulting  in  efficient  direct  numerical  simulation  capabilities 
for  nighly  unsteady  3~D  viscous  flows.  Memory  management  considerations  ar 
currently  being  investigated  using  a  two  pronged  approach. 


Direct,  inversion  methodologies  trade  memory  utilization  (a  relatively 


inexpensive  resource)  for  CPU  efficiency  (still  a  scarse  resource).  For 
fine  3-D  grids  (33x33x33)  to  (65x65x65)  the  coefficient  matrices  employed  to 
gain  the  efficiency  of  the  direct  inversion  procedures  exceed  the  in-core 
memory  conveniently  obtainable  on  the  CYBER  205  at  NASA/Langley.  A 
two-pronged  approach  is  being  pursued  to  circumvent  this  difficulty.  First, 
the  strong  numerical  stability  of  the  direct  solution  procedures, 
demonstrated  in  Osswald  et  al.  [1987],  permits  the  storage  of  these 
coefficient  matrices  in  half-precision  (32  bit)  words  rather  than  the  full- 
precision  (64  bit)  words  currently  employed.  This  effectively  doubles  the 
in-core  storage  conveniently  obtainable  on  the  CYBER  205  system  for  holding 
the  coefficient  matrices  without  further  recoding.  In  addition,  use  of  the 
concurrent  I/O  capability  of  the  CYBER  205  system  is  being  investigated  with 
the  concept  of  implementing  program  controlled  concurrent  paging  through 
this  data  base  simultaneously  with  the  solution  calculations.  This  could 
greatly  reduce  the  in-core  storage  requirement  of  the  direct  solution 
techniques  placing  the  burden  instead  upon  the  concurrent  I/O  processor. 
Finally,  access  to  the  CRAY-2  facility  at  NASA/Ames  is  being  pursued  where 
significantly  more  memory  is  available  to  the  single  user. 

Results  of  this  memory  management  exercise  as  well  as  fine  grid  flow 
results  for  the  driven  cavity  and  3~D  backstep  channel  geometries  has  been 
submitted  for  presentation  at  the  ICNMFD  '83  conference  at  Williamsburg, 
Virginia.  The  flow  inside  a  backstep  channel  with  finite  spanwise  extent 
will  be  examined  next.  For  this  geometry,  careful  1 aser-Doppl er  velocimetry 
results  are  available  for  the  entire  range  of  flow,  i.e.,  for  laminar, 
transition  and  turbulent  flows. 
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3.  RESEARCH  FORECAST 


As  described  in  the  preceding  section,  both  areas  of  research  pursued 
have  progressed  well  to  meet  the  objectives  of  this  study.  During  the 
second  year  of  the  grant  period,  i.e.  October  15,  1987  -  October  14,  1988, 
effort  will  continue  as  follows: 

a.  The  active-flow  reattachment-control  study  will  be  completed  by  carrying 
out  a  parametric  study  of  the  backstep  channel  flow  with  control  flap  to 
achieve  optimum  control.  This  effort  will  be  coordinated  with  that  of 
Nagib  at  the  Illinois  Institute  of  Technology.  Various  scenarios  of 
flap  geometry,  location  and  motion  will  be  examined  in  order  to 
determine  the  optimum  configuration  of  this  flow-control  mechanism. 
Control  of  unsteady  separated  flow  will  also  be  examined  via  the 
elliptic  cylinder  configuration  with  an  unsteady  free  stream. 

b.  Concurrently,  the  3~D  incompressible  unsteady  flow  analysis  and  solution 
procedure  will  be  applied  to  the  flow  in  a  rectangular  cross-section 
duct  with  a  backward  facing  step  and  verified  by  comparing  the 
predictions  with  the  corresponding  1 aser-Doppl er  velocimetry  (LDV)  data 
available  from  the  University  of  Karlsruhe  and  Stanford  University  in 
order  to  obtain  insight  into  the  mechanism  of  shear-layer  instability, 
transition  and  the  strange-attractor  theory  of  turbulence. 

Effort  will  also  be  directed  towards  gaining  better  familiarity  with 
effective  usage  of  the  IRIS  workstation.  This  equipment  is  especially 
useful  in  the  present  unsteady  flow  research,  particularly ,  the  3_D  work. 
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APPENDIX 

EFFICIENT  GENERATION  OF  CONFORMAL  COORDINATES  FOR  ARBITRARY  GEOMETRIES 


Outline  of  Conformal  Mapping  Employed 


Following  the  work  of  Davis  [1983J,  the  Schwar z-Christof f el  conformal 
transformation  given  by  the  equation 


.  n  a.  /it 

rr  M  ^  1 


(A.  1  ) 


is  employed  to  map  an  n-sided  polygon  in  the  z- plane  onto  the  upper  half  of 
the  5- plane  as  shown  in  Fig.  A.1.  In  Eq .  ( A . 1 )  ,  the  a.  ' s  are  known  from  the 

given  geometry,  while  M  and  the  a. ' s  are  unknown. 

Since  it  is  more  appropriate  to  transform  an  arbitrary  channel  to  a 
straight  channel  rather  than  to  an  upper-half  plane,  a  further 
transformation  is  introduced  which  conformally  maps  the  upper-half  5- plane 
to  a  straight  channel  in  the  T- plane.  Clearly,  this  latter  transformation 
is  simply  a  special  eas^  of  the  Schwar z-Chri stoff el  transformation  [Eq. 
(A.1)]  and  is  gi ven  as 
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so  that 
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Hence,  the  overall  transf ormat ion  is  obtained  as 


dT  ~  715  dc 


C=exp  (  tt T ) 


with  dz/d£  given  by  Eq .  (A.1).  Further  details  associated  with  the  actual 
determination  of  this  transformation  are  presently  being  documented  in  a 
paper  (U.  Ghia,  Zuo  and  K.  Ghia  [1988])  to  be  submitted  for  publication  in  a 
technical  journal.  This  work  also  constitute  the  basis  of  a  Master's  degree 
thesis  (Zuo  [1987  ])  . 

Figure  A. 2  shows  typical  grids  obtained,  using  this  analysis,  for  the 
backstep  channel  with  an  oscillating  flap.  Two  representati ve  positions  of 
the  flap  are  shown.  The  grid  clustering  occurring  on  the  left  of  the  flap 
is  unnecessary  for  these  positions  of  the  flap.  However,  it  is  retained 
because  the  gr i d-cluster i ng  transf ormat ions  are  time-invariant  and  the 
clustering  shown  will  be  needed  when  the  flap  rotates  to  occupy  a  position 
in  the  second  quadrant  with  reference  to  the  axis  of  the  flap  motion.  This 
is  a  small  price  to  pay  for  the  important  advantage  that  the  stream-function 
equation  retains  time-invariant  coefficients  and  can  continue  to  be  solved 
efficiently  by  the  BGE  direct  method. 

A . 2  Enhancement  of  Efficiency  for  Evaluation  of  Transformation 

Since  the  transformation  needs  to  be  re-evaluated  for  every  incremental 
change  in  orientation  of  the  flap,  considerable  effort  was  focused  on 
enhancing  the  efficiency  of  determination  of  the  transformation.  The 
results  of  this  effort  are  summarized  in  Table  A.1.  Prior  to  commenting  on 
these  results,  it  is  important  to  note  the  following.  The  Schwarz- 
Christoffel  transf ormat i on  directly  provides  the  complex  derivative  dz/dt  or 
dz/dT  and,  hence,  the  scale  factor  h  of  this  conformal  mapping.  The  flow 
equations  contain  only  h  and  its  derivatives.  Therefore,  integration  of 
dz/dc  itself  is  needed  only  when  it  is  desired  to  view  the  grid.  This  may 
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not  be  necessary  at  each  and  every  time  step  of  the  flow  simulation;  only  h 
needs  to  be  evaluated  at  every  computational  time  level. 

For  channels,  Eq .  (A. 3),  together  with  Eq .  (A.11),  possesses  a  simple 
asymptotic  form  as  Re(c)  -*■  ±  ®.  These  equations  approach  this  asymptotic 
form  in  an  exponential  manner,  so  that  the  simplified  form  constitutes  an 
accurate  representation  of  the  actual  transformation  beyond  relatively  small 
distances  away  from  the  last  deviation  from  the  straight-channel 
configuration.  As  shown  in  Table  A.1,  implementation  of  this  asymptotic 
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form  resulted  in  66-80  percent  saving  of  CPU  time  for  the  computation  of  h  , 
while  the  CPU  time  saved  in  the  grid  calculation  was  almost  80-95  percent 
when  the  grid  is  obtained  by  integrating  dz/dT  in  the  T-plane,  i.e.,  in  the 
straight-channel  plane.  If  the  integration  is  done  in  the  computational 
plane  of  the  clustered  coordinates  (S,N),  the  CPU  time  saved  is  not  nearly 
as  large,  but  the  grid  clustering  obtained  in  the  physical  plane  is  better. 

In  Table  A.1,  item  B  refers  to  the  actual  coding  employed  for 
evaluating  the  mul ti -exponent i al  terms  appearing  in  Eq.  (A.1).  After 
examining  the  manner  in  which  such  terms  are  actually  evaluated  in  the 
computer,  the  coding  was  modified  to  reflect  this  manner  of  evaluation. 

This  resulted  in  an  additional  50  percent  saving  in  the  CPU  time  required  to 
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evaluate  h“ .  Hence,  the  implementation  of  both  items  A  and  B  resulted  in  a 
total  CPU  time  saving  of  almost  90  percent.  It  is  believed  that  this  level 
of  efficiency  is  acceptable  for  proceeding  on  with  the  flow-control  study. 

A  paper  based  on  this  efficient  grid-generation  study  is  presently 
under  preparation  by  U.  Ghia  et  al .  [1988].  This  study  should  also  be  very 
useful  in  various  related  areas  such  as  those  occurring  in  problems 
involving  fluid-structure  interactions  and  free  surfaces,  metal- formi ng 
problems  as  well  as  in  bio-fluid  mechanics  involving  flow  through  passages 
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with  flexible  walls.  All  of  these  problems  encounter  temporally  deforming 
boundaries  and  may  often  require  the  consideration  of  more  than  one 
discipline  of  study.  For  example,  in  fluid-structure  interaction  problems, 
the  fluid  dynamics  analysis  provides  the  time- varying  force  distribution, 
while  the  structural  analysis  provides  the  boundary  deformations 
(vibrations)  occurring  under  the  action  of  these  unsteady  forces.  The  high 
efficiency  of  the  grid-generation  method,  together  with  the  present 
formulation  of  the  stream- fun ction  equation  with  constant  coefficients  even 
in  the  presence  of  temporally  deforming  geometry,  paves  the  way  for 
effectively  examining  the  oseillating-flap  flow-control  mechanism. 
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TABLE  A-l 


ENHANCEMENT  OF  EFFICIENCY 


BY  IMPLEMENTING 
ITEM  'A' 

ITEM  ’B’ 

ITEMS  ’A'+'B’ 


PERCENTAGE  OF  COMPUTATIONAL  TIME  SAVED 


h  CALCULATION 


GRID  CALCULATION 


66  -  80 
50 

83  -  90 


80  -  95*  (66  -  80) 
0 

80  -  95*  (66  -  80) 


+ 


+ 


'A'  :  ASYMPTOTIC  FORM  FOR  CHANNEL  MAPPING. 

’B'  :  TIME-SAVING  TECHNIQUE  FOR  EVALUATING 

MULTI-EXPONENTIAL  TERMS. 

NOTE:  CURVED  BOUNDARIES  ARE  ALSO  TREATED  WITH 
COMPARABLE  EFFICIENCY. 


*  INTEGRATION  DONE  IN  CHANNEL  PLANE  (T). 

+  INTEGRATION  DONE  IN  COMPUTATIONAL  PLANE  (S,N) 
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FIG.  A-l  CONFORMAL  MAPPING  OF  ARBITRARY  CHANNEL 
TO  STRAIGHT  CHANNEL. 


USTERBD  CONFORMAL  BOUNDARY-ALIGNED  GRID  FOR  BACKSTEP  CHANNEL 
TH  OSCILLATING  FLAP  IN  TWO  DIFFERENT  POSITIONS. 


